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Abstract. The NA49 experiment at CERN SPS has acquired a huge data set of 
Pb+Pb events over a broad range of energy and centrality during the last several years. 
This high statistics data set, coupled with a state-of-the-art analysis technique, allows 
for the first model-independent extraction and energy scan of 3D emission sources 
for pion pairs at SPS energies. These 3D pion emission sources provide new insights 
into the nature of a long-range source previously reported by PHENIX at RHIC. The 
new results indicate that the pion source displays significant non-Gaussian tails in the 
longitudinal direction at 40 and 158 AGeV and in the outward direction at 158 AGeV. 



PACS numbers: 25.75.-q, 25.75. Gz, 25.70.Pq 
1. Introduction 

A deconfined phase of nuclear matter is expected to be formed at the high energy 
densities created in relativistic heavy ion collisions [lj. It is widely believed that 
important signatures of such a phase are reflected in the space-time extent and shape 
of particle emission source functions. 

Recently, 1-Dimensional source imaging techniques [4J [5] have revealed a non-trivial 
long range structure in the two-pion emission source at RHIC (2j [3] . The origins of this 
structure are still unclear. The presence/absence of such a structure in the pion emission 
sources in heavy ion collisions at intermediate SPS energies could yield important 
information which could help resolve the structure's origins. The NA49 Collaboration 
has carried out Pb+Pb collisions over a wide range of bombarding energies at the CERN 
SPS during the last decade [6|. Such a rich data set provides a unique opportunity to 
search for this long range structure at the SPS and study its evolution with beam energy 
with a view to unraveling its nature. 

In this paper, the 3-Dimensional emission source images for pions produced in 
central Pb+Pb collisions over the incident energy range 40 and 158 AGeV are presented. 
The results are discussed in the context of a Gaussian shape assumption. 
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2. Experimental Setup and Data Analysis 

The data presented here were taken by the NA49 Collaboration during the years 1996- 
2002. The incident beams of 40 and 158 AGeV were provided by the CERN SPS 
accelerator. The NA49 Large Acceptance Hadron Detector [7] achieves large acceptance 
precision tracking (Sp/p 2 ~ (0.3 — 7).10~ 4 (GeV/c)~l) and particle identification using 
time projection chambers. Charged particles are detected by tracks left in the TPC and 
identified by the energy deposited in the TPC gas. Mid-rapidity particle identification 
is further enhanced by a time-of- flight wall (resolution 60ps). Event centrality is 
determined by a forward calorimeter which measures the energy of spectator matter. 

3D correlation functions , C(q), were obtained as the ratio of pair to uncorrelated 
reference distributions in relative momentum q for it~it~ pairs. Here, q = teizEal is half 
of the relative momentum between the two particles in the Pair Center-of-Mass System 
(PCMS). The pair distribution was obtained using pairs of particles from the same 
event and the uncorrelated distribution was obtained by pairing particles from different 
events. Track merging and splitting effects were removed by appropriate cuts on both 
the pair and uncorrelated distributions. Momentum resolution effects were negligible. 

In the cartesian harmonic decomposition technique [SJ 12], the 3D correlation 
function is expressed as 

C(q) - 1 = i?(q) = E E RL-aMX-a^) (1) 

I Oil— Oil 

where I — 0, 1, 2, . . ., aij = x, y or z, A l ai (flq) are cartesian harmonic basis elements 
(Qq is solid angle in q space) and R ai a[ (q) are cartesian correlation moments given by 

R l ai ... ai (i) = { ^r^ J ^<..*n(&d m (2) 

The coordinate axes are oriented so that z is parallel to the beam (long) direction, x 
points in the direction of the total momentum of the pair in the Locally Co-Moving 
System (LCMS) (out) and y is perpendicular to the other two axes (side). 

The correlation moments, for each order I, are calculated from the measured 3D 
correlation function using equation (j2j). Each independent correlation moment is then 
imaged using the ID Source Imaging code of Brown and Danielewicz [31 IS] to obtain 
the corresponding source moment for each order I. Bose-Einstein symmetrisation and 
Coulomb interaction (the sources of the observed correlations) are contained in the 
source imaging code. Thereafter, the total source function is constructed by combining 
the source moments for each I as in equation (jHJ) 

%) = E E i... ai M<... a ,(n r ) (3) 

l ai—at 

3. Results 

Figure [T] shows the I = (R°) and I = 2 (Rl 2 an d Ry 2 ) moments for mid-rapidity 
(\Vl — Vo\ < 0.35, where jjl is particle laboratory rapidity and y is CM rapidity), lowp^ 
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Figure 1. 1 = and I = 2 moments for mid-rapidity low px tt~-k~ pairs 
from 158 AGeV central Pb+Pb collisions as a function q. Data are shown as 
solid circles while the squares represent the result of a simultaneous fit of the 
moments with an ellipsoid shape ((a)-(c)) and a 2-source model ((d)-(f)). 



(0 < pt < 70MeV/c) 7r~7r~ pairs from 158 AGeV central (7%) Pb+Pb collisions, as a 
function of the relative momentum q in the PCMS frame. Here, R^. 2 is shorthand for R^ x 
etc. Moments for other l's are either zero or negligible. The Lorentz transformation of q 
from the laboratory frame to the PCMS is done by transforming to the LCMS along the 
beam axis and then transforming to the PCMS along the pair transverse momentum. 

The data are represented in Fig. [T] by solid circles, while squares represent the 
results of fits to the independent moments with model sources, an ellipsoid shape 
(3D Gaussian) in panels (a)-(c) and a 2-source shape in panels (d)-(f). The last 
function arises by assuming a linear combination of two Gaussians for the single particle 
distribution. The best-fit parameters for the two shapes are listed in the lower panels. 
The 2-source shape yields a much better representation of the data. 

Given that the only significant imaged moments are found for I = and 2 
multipolarities, the net imaged source function in the x, y and z directions is simply the 
sum of the ID source S° and the corresponding I = 2 moment Sf t where i= x, y or z. 
In panels (a)-(c), Figure [2] compares the net imaged source (squares) for mid-rapidity, 
low pt pion pairs from central Pb+Pb collisions at 40 AGeV (i.e ^/s = 6.4AGeV) to 
the best-fit parametrized functions, Gaussian (triangles) and 2-source model (circles). 
The source image and the 2-source model agree very well in all 3 directions and both 
disagree with the Gaussian fit in the z direction. 

At 158 AGeV (i.e = 17.3AGeV), represented in panels (d)-(f) of Fig. [21 the 
findings are similar, but in addition to the z direction, the source function starts to 
develop a non-Gaussian tail in the x direction. 
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Figure 2. Source function in x (out), top panels, y (side), middle panels, and z 
(long) directions for low px mid-rapidity 7r - 7r - pairs from central Pb+Pb events 
at yfs = 8.9AGeV (incident beam energy of 40 AGeV) and yfs = 17.3AGeV 
(incident beam energy of 158 AGeV). The imaged, Gaussian and 2-source 
functions are represented, respectively, by squares, circles and triangles. 



4. Discussions 

The ratio of the RMS radii of the source functions in the x and y directions is 1.3±0.1 at 
40 AGeV and 1.2±0.1 at 158 AGeV. This deviation from unity, evident visually, points 
to a finite pion emission time. 

Moreover, the RMS pair separation in the z direction, from figure [2]^ c) and (f), is 
llfm at 40 AGeV and 12fm at 158 AGeV. These dimensions are much smaller than the 
Lorentz-contracted nuclear diameters of 3fm at 40 AGeV and 1.5fm at 158 AGeV and 
in fact give the RMS pair separation due to the longidutional spread of nuclear matter 
created by the passage of the two nuclei. Since the latter are moving with almost the 
speed of light, one can infer the lower bound formation time of the created nuclear 
matter to be 8fm/c at 40 AGeV and lOfm/c at 158 AGeV. 
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